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Abstract 

We comparatively analyze the rare Aj, —>• channel in standard model, su¬ 

persymmetry and Randall-Sundrum model with custodial protection (RS C ). Using 
the parametrization of the matrix elements entering the low energy effective Hamilto¬ 
nian in terms of form factors, we calculate the corresponding differential decay width 
and lepton forward-backward asymmetry in these models. We compare the results 
obtained with the most recent data from LHCb as well as lattice QCD results on the 
considered quantities. It is obtained that the standard model, with the form factors 
calculated in light-cone QCD sum rules, can not reproduce some experimental data 
on the physical quantities under consideration but the supersymmetry can do it. The 
RS C model predictions are roughly the same as the standard model and there are no 
considerable differences between the predictions of these two models. In the case of 
differential decay rate, the data in the range 4 GeV 2 /c 4 < q 2 < 6 GeV 2 /c 4 can not 
be described by any of the considered models. 

PACS number(s): 12.60.-i, 12.60.Jv, 13.30.-a, 13.30.Ce, 14.20.Mr 


*e-mail: kazizi@dogus.edu.tr 
i e-mail: tugba.olgun@okan.edu.tr 

1 e-mail: zeynep.tavukoglu@okan.edu.tr 



1 Introduction 


The ATLAS and CMS Collaborations at CERN have independently reported their discovery 
of the Higgs boson with a mass of about 125 GeV using the samples of proton-proton 
collision data collected in 2011 and 2012, commonly referred to as the first LHC run [1-3]. 
Recently, a measurement of the Higgs boson mass based on the combined data samples of the 
ATLAS and CMS experiments has been presented as rrin = 125.09 ± 0.21 (stat) ±0.11(s|/st) 
GeV in Refs [3-7]. At the same time, all LHC searches for signals of new physics above the 
TeV scale have given negative results. However, the LHC constraints on new physics effects 
can help theoreticians in the course of searching for these new effects and answering the 
questions that the standard model (SM) has not answered yet. We hope that the upcoming 
LHC run can bring unexpected surprises to observe signals of new physics in the experiment 
[ 8 ], 

Although the SM could be valid up to some arbitrary high scale, new scenarios should 
exist because we are lacking a proper understanding of some important issues like origin 
of the matter, matter-antimatter asymmetry, dark matter and dark energy etc. [9]. In 
the baryonic sector, the loop-induced flavor changing neutral current (FCNC) decay of the 
Af, —>■ A i + £~ with l = e,n,r, which is described by the b —> s£ + £~ transition at quark 
level, is one of the important rare processes that can help us in the course of indirectly 
searching for new physics effects [10]. Recently, the differential branching fraction of the 
A® —> A/j, + /j,~ decay channel has been measured as a function of the square of the di-muon 
invariant mass (g 2 ), corresponding to an integrated luminosity of 3.0 /5 _1 using proton- 
proton collision data collected by the LHCb experiment [11]. The measured result at 15 
GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 region for the differential branching fraction is dBr(A° —y 
A fi + ij,~)/dq 2 = (1.18 ± 0.03 ± 0.27) x 10 - ' GeV 2 /c 4 . The LHCb Collaboration has 

also reported the measurement on the forward-backward asymmetries of this transition at 
the // channel. The measured result at the 15 GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 region for the 
lepton forward-backward asymmetry is A BB = —0.05 ± 0.09 (stat) ± 0.03 (syst) [11]. In the 
literature, there are a lot of studies on this decay channel via different approaches (for some 
recent studies see for instance Refs. [12-17]). 

In the present work, we calculate the differential decay rate and lepton forward-backward 
asymmetry related to the FCNC A h —> A£ + £~ transition for all leptons in the SM, super- 
symmetry (SUSY) and Randall-Sundrum scenario with custodial protection (RS C ). We 
compare the results with the experimental data provided by LHCb [11] as well as the ex¬ 
isting lattice QCD predictions [18]. Comparison of the LHCb results with the lattice QCD 
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predictions shows that there are some deviations of data from the SM predictions. Such 
deviations can be attributed to the new physics effects that can contribute to such loop level 
processes. In this connection, we comparatively analyze the A& —y A £ + t~ decay channel in 
SM and some new physics scenarios. In the calculations, we use the form factors calculated 
via the light-cone QCD sum rules in [19]. Hence, to get ride of any misleading, we will 
use SMLCSR instead of SM referring to the results that are obtained via using the from 
factors predicted by the light-cone sum rules in [19] when we speak about the predictions of 
different models. Note that there are many studies devoted to the calculations of the form 
factors defining the transition under consideration via different approaches (se for instance 
[ 20 , 21 ]), but our aim here is to use those form factors that are obtained in the full theory 
of QCD in [19] without any approximation. 

The outline of the paper is as follows. In the next section, we introduce a detailed 
discussion of the effective Hamiltonian responsible for the semileptonic A;, —» A£ + £~ decay 
channel and Wilson coefficients in SM, RS C and SUSY models. In this section, we also 
present a basic introduction of the RS C scenario. In section 3, we calculate the differential 
decay rate and lepton forward-backward asymmetry at different scenarios and compare the 
predictions of different models. 


2 The semileptonic —> A£ + £ transition in SM, 

SUSY and RS C models 


2.1 The effective Hamiltonian and Wilson Coefficients 


At quark level, the FCNC transition of A & —» A£ + £ is governed by the b —> s£ + £ transition 
whose effective Hamiltonian in the SM can be written as 


-1ye// 
n SM 


G /Vr ,111 1 / b 1 / 


ts 


2vA 


7T 


Cg ff ’ SM S ^u(l - ^)b£^£ + Cf 0 M s 7 M (l - ^)b£Yld 


- 2m h Cj ff ’ SM —sicr^ w q u (l + 75 )b i^i 

qz 


( 2 . 1 ) 


where and V t * are elements of the Cabibbo-Kobayashi-Maskawa (CKM) mixing ma¬ 
trix, a em is the fine structure constant at Z mass scale, Gp is the Fermi weak coupling 
constant, q 2 is the transferred momentum squared; and the Cg ^’ SM , Cf 0 M and C e ^' SM 
are the Wilson coefficients representing different interactions. The explicit expressions of 
the Wilson coefficients entered to the above Hamiltonian are given in the following. The 
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Wilson coefficient Cg^’ SM which is a function of s' — with q 2 lies in the allowed region 

o 

4 mf < q 2 < (m\ b — m a ) 2 is given by [22, 23] 

C e 9 ff ' SM (s') = C? DR ri(s') + h(z, s') (3Ci + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 


■“7^(1 • s') (4C',3 + 4C*4 + 3C*5 + Cg) 

—-/r(0, s') (C 3 + 3 C 4 ) + - (3C*3 + C 4 + 3C*5 + Cg) , 


( 2 , 2 ) 


where the C 9 DR in the naive dimensional regularization (NDR) scheme is expressed as 

ySM 


qNDR _ pNDR 


sin 2 e 


4 Z SM + P e E s 


'w 


(2.3) 


The last term in the right hand side is neglected due to smallness of the order of P E . Here 
P 0 NDR = 2.60 ±0.25, Y SM = 0.98, Z SM = 0.679 and sin 2 6 W = 0.23 [22-24], The parameter 
r)(s') in Eq.(2.2) is given as 


7 7 (s') = 1 + ^M u( s>) 


with 


7T 


2 4 2 5 -t- 4 s' 

^ = '9’ r ^3 Li2(l ' ) “3 lllS ' ln(1 “ l ' ) ^yTT2F) ln(1_l ' ) “ 


(2.4) 


2s'(l±s')(l -2s') 


In s' + 


5 ± 9s' - 6 s' 


i/2 


and 


3(1 - s') 2 (l ± 2s') 111 ° ^ 6(1 - s')(l ± 2s') ’ 

a s (m z ) 


(2.5) 


a s (x) = 


aAmz) ln (m z ~ 


1 “ 00 2 n 

Here a s (rriz) = 0.118 and f3o = yr. The function h(y,s') in Eq.(2.2) is also defined by 


( 2 . 6 ) 


, / -/\ 8 m b 8 8 4 

ft «, s = — In-In y H- —x 

Ky ' J 9 n h 9 y 27 9 


(2.7) 


In 


/T^x+i 


--{2 + x )\l~ x \ 1/2 
y 1 2 arctan 


l-x-l 

1 


/x—1 ' 


• in j , for x = < 1 

for x = > 1 , 


( 2 . 8 ) 


where y = 1 or y = z = — and, 

& & mj, 7 


h( 0, s') = 


27 


8 , rn b 

- In — 

9 y b 


4 1 v , 4 • 
9 lnS W 


(2,9) 
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In Eq.(2.2), the remaining coefficients are given by [24] 


8 

Cj = £ 

2=1 

where the kji are given as 


hi = ( 

0 , 

0 , 

1 

2’ 

1 

2’ 

0. 

hi = ( 

0 , 

0 , 

1 

2’ 

1 

2’ 

0. 

k 8 i ( 

0 , 

0 , 

1 

14’ 

1 

6’ 

0.0510 

hi = ( 

0 , 

0 , 

1 

14’ 

1 

6’ 

0.0984 

hi = ( 

0 , 

0 , 

0, 

0, 

-0.0397 

hi = ( 

0 , 

0 , 

0, 

0, 

0.0335 


(.7 = 1,- 

.6), 


(2.10) 

0 , 

0 , 

0 ) , 


0 , 

0 , 

0 ) , 


-0.1403, 

-0.0113, 

0.0054 ) , 

(2.11) 

0.1214, 

0.0156, 

0.0026 ) , 

0.0117, 

-0.0025, 

0.0304 ) , 


0.0239, 

-0.0462, 

-0.0112 ) . 



The explicit expression for the Wilson coefficient Cf 0 M is given as 

ySM 


r<SM _ 
L '10 


sir 6 W 


( 2 . 12 ) 


Finally, the Wilson coefficient C 7 ^ ,SM in the leading log approximation is defined by [22-25] 

8 8 

= r)^C 7 {n w ) + - (7735 - r)^j C 8 (nw) + C 2 (nw) ^ hr) ai , 


2=1 


(2.13) 


where 


V = 


Qq(hw) 


(2.14) 


and 


CV(hrn) — ~ 2 D'o SM ( x t) , C 8 {^w) — — ~E' 0 SM {x t ) , C 2 {^w) — 1 • 


(2.15) 


The functions D' 0 SM (xt) and E' 0 SM (xt) with Xt = —h are given by 


D' 0 SM (x t ) = 


(8x% + bx 2 — 7x t ) x 2 (2 — 3x t ) 


+ 


12(1 — x t y 2(1 — x t ) 


In x t , 


(2.16) 


and 


K SM (x t ) = 


x t (x 2 - 5x t - 2) 


3x 2 

+ —-rr hi x t . 


4(1 ~x t ) 3 2(1-ay) 4 


(2.17) 


4 



The coefficients h t and a* inside the C^' SM are also given by [22,23] 

hi = ( 2.2996, -1.0880, -f, --k, -0.6494, -0.0380, -0.0186, -0.0057 ), 

I 

J., -||, 0.4086, -0.4230, -0.8994, 0.1456 ). 


a,- = 


14 

23’ 


(2.18) 


One of the most important new physics scenarios is SUSY. The different SUSY models 
involve the SUSY I, SUSY II, SUSY III and SUSY SO(10) scenarios according to the values 
of tanfd and an extra parameter // with dimension of mass [26-29]. In SUSY I, the Wilson 
coefficient C * 7 changes its sign, the /i takes a negative value and the contributions of the 
neutral Higgs bosons (NHBs) have been disregarded. In the SUSY II model, the value of 
the tan[3 is large and the masses of the superparticles are relatively small. In SUSY III, 
the tanfd takes a large value and the masses of the superparticles are relatively large. In 
SUSY SO(10), the contributions of the NHBs are taken into account. The supersymmetric 
effective Hamiltonian in terms of the new operators coming from the NHBs exchanges 
diagrams and the corresponding Wilson coefficients is written as 


nj e ff 

rt SUSY 


GFCUemVtbVt 


ts 


2 y/ 2 - 


TT 


ci ff ’ SUSY s llx ( i - l5 )bhn + c' eff ’ SUSY s^( i + 75 )bhn 


+ c 


SUSY 
10 


S7#i ( 1 ^ 75)6 + C[ 0 SUSY s'y li (l + 75)6 £ 7 ^ 75 ^ 


- 2 m b C? ff ' SUSy —sicr^q u (l + 75 )b t^t - 2 m b C' 7 eff ’ 
+ C™ SY s( 1 + l5 )bU + C' q * usy s( 1 - -K)bU 

+ C s q u 2 sy s (1 + 75)^75^ + C ' q s 2 usy s ( 1 - 75)6*75* 


SUSY 


-sia^q v {l - 75 ) 6 ^ 7 ^ 


(2.19) 


sc/5y 


r eff,SUSY r t eff,SUSY r SUSY rv SUSY r eff,SUSY w eff,SUSY r SUSY rv S 

wnere o 9 , o 9 , o 10 , o 10 , o 7 , o 7 , Oq 1 , Oq 1 

Cq^ s ^ and Cq^ usy are the new Wilson coefficients in the different SUSY models. The new 
Wilson coefficients, Cq^ US ^ and Cq 2 SUS ' 1 come from NHBs exchanging [29]. The primed 
Wilson coefficients only appear in SUSY SO(10) model. The values of Wilson coefficients 
in different supersymmetric models are presented in table 1 [27-30]. 

The last new physics scenario which we consider in this work is the Randall-Sundrum 
scenario proposed to solve the gauge hierarchy and the flavor problems in 1999 [31,32]. It 
is a successful model, featuring one compact extra dimension with non-factorizable anti- 
de Sitter ( AdS 5 ) space-time [33]. This model describes the five-dimensional space-time 
manifold with coordinates (x; y) and metric 


ds 2 = e 2kVr q flu dx IJ 'dx l/ — dy 2 , 
Vu* = diag(+ 1 ,- 1 ,- 1 ,- 1 ) . 


( 2 . 20 ) 
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Coefficient 

SUSY I 

SUSY II 

SUSY III 

SUSY SO(10)(Y 0 = -1000) 

c eff,SUSY 

+0.376 

+0.376 

-0.376 

-0.219 

c eff,SUSY 

4.767 

4.767 

4.767 

4.275 

rsusY 

-3.735 

-3.735 

-3.735 

-4.732 

csusy 

U Q i 

0 

6.5(16.5) 

1.2(4.5) 

0.106+ 0i(l.775+ 0.002i) 

r<SUSY 

0 

—6.5(—16.5) 

—1.2(—4.5) 

-0.107 + 0i(—1.797 - 0.002*) 

(j! eff,SUSY 

0 

0 

0 

0.039 + 0.038i 

r t eff,SUSY 
°9 

0 

0 

0 

0.011 + 0.072i 

r" susy 
°10 

0 

0 

0 

-0.075 - 0.67i 

rv SUSY 

L 'Q i 

0 

0 

0 

-0.247 + 0.242i(—4.148 + 4.074i) 

r» susy 
L 'Q2 

0 

0 

0 

-0.25 + 0.246i(—4.202 + 4.128*) 


Table 1: The Wilson coefficients in different SUSY models [27-30]. The values inside the 
parentheses are for the r lepton. 


The scale parameter k is defined as k ~ 0(M P i anck ). We choose it as k — 10 19 GeV. The 
fifth coordinate y varies in a range between two branes 0 and L. y = 0 and y = L correspond 
to the so-called UV brane and IR brane, respectively. The simplest RS model with only the 
SM gauge group in the bulk has many important problems with the clectroweak precision 
parameters [34]. In the present work, we consider the RS model with an enlarged custodial 
protection based on SU( 3) c x SU{2)l x SU(2)r x 17(1) x x Plr, where Plr interchanges 
the two SU( 2) groups and is responsible for the protection of the Zb^L vertex (for more 
information on the model see [33-41]). 

The effective Hamiltonian for the b —> s£ + £~ transition in the RS C model is given as 


n 


eff 

RS C 


G FOt-emVtbVts 
2y/2lX 


ci f f'RScs llx (i - lb )bhn + c' g eff ' RSc s lfl (i + lb )bhn 


+ C^s^l - lb )blnd + c[ 0 RS ‘s^(l + ^btYld 

- 2m b C 7 ff ' RSc ^sia^q u (l + y 5 )bl^i 


- 2m b C' 7 ef k RS -}^ s i a q ^[i _ y 5 )fe7y^7 , 

qz 


( 2 . 21 ) 


where the new Wilson coefficients are modified considering the new interactions. The new 
coefficients in terms of the SM coefficients are written as [33-41] 


fj(i)RS c 


C 


(. l)SM 


+ a c. 


(b 


i = 7,9,10 , 


( 2 . 22 ) 
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where 


and 


with 


and 


AC 9 

ACg 

ACio 


AK 


siri 2 (6 w ] 

A Y[ 

sin 2 {6 W/ 

AK 

sin 2 [9 w 


-4A Z, 


- 4AZ' 


(2.23) 


An 

An' 

az s 

AZ' S 


ACJo = 

An 

sin 2 (9 w ) 

1 V 

Af(A') - A«(.A) 

v tb v t * s ^ 

^Mx9sm 

1 V 

A“(.A) - A«(.A) 

v*vz ^ 

4 Mx9sm 

1 \n 

K(x) Ai 

V tb V t * s ^8M 2 x g 2 SM sin 2 (9 w ) 

1 V' 

Ag(X) 


^l(X) , 
Ar(A) , 
A$f(A) , 


Xr(X) . 


v tb v t * s ^ m 2 x gl M sin 2 {e u 


(2.24) 


(2.25) 


(2.26) 


In the above equations, X = Z, Z#, Z' and g| M = ^§ 2 tt sin?(e ) an< ^ ^ ie Weinberg 
angle. The functions inside An, A Y', AZ s and A Z' s are given in [33-41], 

In the case of A C%\ ACj\g b ) = 0.429AC'y) (Mata;) + 0.128 AC%\Mkk) is used where 
the following three contributions are included [35]: 


(ACV)i 

(yACj^2 


(A C 7 ) 3 


iQ u r [A + 2m 2 F {A' + B')]\v{Y u {Y u )^Y d V R 


F=u,c,t 
8 
3* 


23 


—iQd r~(gt D y [n + A + B + 4m 2 F (I' 0 + A' + B')\ 

F=d,s,b 

v\n L Y d n R v R } 


23 

iQdr^{g A s D ) 2 rn F [I 0 + A + 

F=d,s,b 

— \v[Y d n R n R v R 

m s L J 23 


v[n L n L Y d v R 


J 23 


(2.27) 
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(AC'), = 
(A C') 2 = 

(A C' 7 ) 3 = 


iQ u r Y ^ + 2 m 2 F {A' + B ')] \p ] R {Y d )^Y u {Y u )^V L ^ j 

F=u,c,t 

-iQd r-(gt D y [lo + A + B + 4m 2 F (I' 0 + A'+ B')] 


F=d,s,b 


v\n R (Yyn L v L 


23 


iQd r-{g A s D ) 2 


: m b 

m s 


Y m F [I 0 + A + B]{ [tf R K R K R (Y d )'V L 

F=d,s,b 

v^Y d yn L n L v L 


23 


23 


(2.28) 


(AC's) 

(AC 8 ) 


i — 


2 — 


ir E [A + 2m 2 F (A' + B ')] \v ] L Y u (Y u )^Y d V R 

F=u,c,t 

■ 9 / 4_D \ 2 77 

-vr-{g s ) - 

8 F=d,s,b 

v\Y d n R (Y i yn L Y d v R 

9 


J 23 


T 3 Y [A + B + 2m 2 F (A' + B')] 


J 23 


(ACs) 3 — ~'i' r ^(gt D ) 2 7 3 ^ [A + B + 2 m 2 F (A' + B')] 

F=d,s,b 


v\n L Y d n R v R 


23 


(2.29) 


(ACg)i = ir Y [A + 2m 2 F (A' + B')}\v^ R {Y d yY u (Y u yV L 

F=u,c,t 


J 23 


(AC') 2 = -!r|( 9 f) 2 — 
8 rrib j m s 


T 3 Y [A + B + 2m 2 F {A! + B')\ 


F=d,s,b 

v ] R (Y d yn L Y d n R (Y d yv L 

9 


23 


(ACg) 3 — -ir-^gF) 2 % Y, [A + B + 2m F (A' + B')\ 

F=d,s,b 


v ] R n R (Y d yn L v L 


23 


where r = 


and % = j J„ L eh/[g(?/)] 2 . For the parameters inside the above equa- 


(2.30) 


^FA s m b “““ L JO 


tions and the related diagrams see [33-41], The Q u and Qd are representing the electric 
charges of the np and down type quarks, respectively. The functions l{f\ A^ and B (l) are 



given as 


m 

w 

A(t ) 
A'(t) 
A(t) 
A'(t) 


(An) 2 Ml K 

i 1 


t - 1 
1 +1 


+ 


ln(t ) 

(^l ) 2 

ln[t) 

(47r) 2 ^ 2t(t — l) 2 (f — l) 3 

Bit) = ' 1 


2 B\t) = 


m = 


(4vr) 2 4 M 2 KK 
i 1 


t — 3 2ln(t) 


(4vr) 2 
i 1 


fi'(t) = 


(4vr) 2 4 M 2 KK 
i 1 


(t-1) 2 (t-l) 3 

t 2 — 5f — 2 /n(t) 

6 t(t - l) 3 ~~ (7 - 1 y 

‘it — 1 2 t 2 ln(t) 


it- 1) 2 (t-1) 3 

5t + 1 2t(2 + t)ln(t ) 


(4vr) 2 4M* * l (t - 1)* 


(t - 1)< 


(2.31) 


with t = ’m 2 F /M‘f <K (for more information see [35]). 

Fitting the parameters to the B —> K*p + p~ channel, the modifications on Wilson 
coefficients in RS C model are found as table 2 [35]. 


A C 7 

A C' 7 

ac 9 

< 

ACio 

A^o 

Values 0.046 

0.05 

0.0023 

0.038 

0.030 

0.50 


Table 2: The values of modifications in Wilson coefficients in RS C model used in the analysis 
[35]. 

2.2 Transition amplitude and matrix elements 

The amplitude of the transition under consideration is obtained by sandwiching the corre¬ 
sponding effective Hamiltonian between the initial and final baryonic states, i.e., 

M A ^ M+e ~ = (A(p A ) I n eff I A 6 (p A J> , (2.32) 

where p\ b and p\ are momenta of the initial and final baryons. To calculate the amplitude, 
we need to know the following matrix elements which are parametrized in terms of twelve 
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form factors in full QCD: 


(A(Pa) I s 7 /,(! - 75 )b | MPA b )) = «a(Pa) 


iMq 2 ) + *<v772 (g 2 ) + 77s (g 2 ) 


~ Inls9i{q l ) ~ icr^ 5 q u g 2 (q 2 ) - 75.93(<? 2 ) 


UA b (PA b ) , 

(2.33) 


(A(pa) | 57^(1 + 75)6 | Ab(p A J) = ua(pa) 


7 / 7 i(g 2 ) + i°^q u f 2 {q 2 ) + g A 7s(g 2 ) 


+ 7m759i(9 2 ) + ivixvi5q v g2{q 2 ) + q^i^g^q 2 ) 


UA b (pA b ) , 

(2.34) 


(A(pa) I Sia lll/ q l '( 1 + 75)6 | A 6 (paJ) = «a(pa) 


lji(q 2 ) + i<y^q v fl (g 2 ) + g A 7 j (g 2 ) 


inisgiiq 2 ) + *<V75 q u 92(q 2 ) + g M 75gJ(g 2 ) 


u A b (PA 6 ) , 

(2.35) 


(A (p A ) I sia, ll/ q l '( 1 - 75)6 | A 6 (p A J) = «a(pa) 


7 /i/i (? 2 ) + ™^q u fl (g 2 ) + g M / 3 T (g 2 ) 


7 ^ 75 Pi (9 2 ) - *<V 75 q u g 2 (q 2 ) - 77 5 gJ(g 2 ) 


u A b (PA b ) , 

(2.36) 


(A(pa) I 5(1 + 75)6 | A 6 (p A J) = —wa(pa) 

m b 


7A(g 2 ) + iq^a^q u Mq 2 ) + g 2 /s(g 2 ) 


- ^7sPi(g 2 ) - iq li Vnul5q u g2{q 2 ) - 9 2 75Ps(g 2 ) 


UA b (pA b ) , 

(2.37) 


and 


(A(pa) I s(l -75)6 | Aj,(p A J) = —m a (pa) 

rn b 


dfi(q 2 ) + iq^v^q 1 ' f 2 (q 2 ) + g 2 /s(g 2 ) 


4ib9i{q 2 ) + 7 M ov 75 <f 92 (g 2 ) + g 2 75 P 3 (g 2 ) 


WpaJ , 

(2.38) 
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where the /'■ T ' > and gf r ' > (i running from 1 to 3) are transition form factors and u Ab and u A 
are spinors of the A h and A baryons, respectively. We will use these form factors from [19] 
that have been calculated using the light-cone QCD sum rules. 

Using the above transition matrix elements in terms of form factors, we find the ampli¬ 
tude of the transition under consideration at different scenarios. In the SM, we find 

= Gf0 £%* K | + B? M L] + ia^f\Al M R + Bl M L) 

+ <f[Ai M R + Bl u L])u^)] (ire) 

+ [s A (PA)(7,[®f M R + £i M L] + ia„„q"[D$ M R + £ 2 SM L] 

+ qH[Dl u R + £i M L])u At (i> A ,)\ (*/W)j • 

(2.39) 

In the case of SUSY we get 

= Gf °!^ K | [uMhMl USr R + Bf SY L] + ,a^[Al USY R + B§ USY L) 

+ q»[A s 3 USY R + B$ USY L])u Ab (p Ab )] {hH) 

+ [u a (pa)(^[Di USY R + Si SUSY L\ + ia^[D s 2 USY R + £ 2 SUSY L] 

+ q»[D s 3 USY R + £$ USY L])u Ab (p A j] 

+ [ua(paMG s i usy R + nf usy L] + iq^^[G! USY R + H S 2 USY L] 

+ q 2 [G° USY R + ni USY L])u Ab (p Ab )\ (U) 

+ [MpaMJC^R + S* USY L] + iq^, v cf[1C s 2 USY R + S s 2 usy L] 

(2.40) 
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and for RS C we obtain 


= Gf °^ K | [s A (p A )( 7 ,[-4f s --R + Bf-L] + iw'lA?‘R + B^L] 

+ <t[A^-R + sf'£])«A*feo] (ire) 

+ [uAtPAKTApP-R + £i S ’L] + ia^lVf’R + £ 2 rs 'L] 

+ <f{Df--R + £ 3 rs 'L])u Ai (paJ] (lr 75«)j , 

(2.41) 


where R = (1 + 7 s )/2 is the right-handed and L = (1 — 7 5 )/2 is the left-handed projectors. 
In the above equations, the calligraphic coefficients are defined at different models as 


A l = /,C S '" + - giC*"- 

- 2m ‘^ fl C 7 ,f+ + 

, >4.2 — *Ai(l —>■ 2), .4-3 — v4i (1 —>• 3) 

Bi = /ic, e//+ + gi c;"- 

- 2m„4 [/rq" + - sfc?"-' 

, S 2 = Si (1 —>• 2), $3 = B\ (1 —>■ 3) 

Di - liC/u - '/iC l(l . 

^2 = Hi (1 2 ) i: 

T>z — T>i (1 —>• 3), 

£1 — /lC+ + 9iC 10 , 

^2 = (1 —^ 2 ) , 

£3 = £\ (1 —» 3), 

Si = — [ /lC + - 9l C 0 j 

m b L V1 V1 J 

, Q 2 — Q\ (1 —>• 2 ) , 

Qi = ^7i (1 —^ 3), 

= ~ f lC Ch + 

,,b b L 

, % = ^ 1 (l-> 2 ), 

Hi = Hi (1 -> 3), 


, K-2 = /Ci (1 —» 2 ), 

/C 3 = /Ci (1 —y 3), 

s ' = ^K + ’ +9 ' c «’] 

, ^2 = (1 —)■ 2 ), 

£3 — £1 (1 — t 3), 


(2.42) 


with 

(ief}+ _ (~< e ff W eff s-ieff- _ s^eff _ W eff 

°9 ~ °9 “T 0 9 i c '9 — L/ 9 °9 ) 
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c e 7 ff+ 

= c e 7 ff + c' 7 eff , 

c e 7 ff - 

= c e 7 ff 

- c 7 eff 

Cto 

= C W 

+ 

c[ o, 

C U 3 

= Ci 0 

- C[ 0 , 

r+ 

L 'Q i 

= c Ql 

+ 

C' Ql > 


= c Ql 

- C' Ql 

r+ 

L Q2 

= Cq 2 

+ 



= Cq 2 

Cq 2 


(2.43) 


3 Physical Observables 

3.1 The differential decay width 

In the present subsection, we would like to calculate the differential decay width for the 
decay channel under consideration. Using the decay amplitude and the transition matrix 
elements in terms of form factors, the supersymmetric differential decay rate as the most 
comprehensive differential decay rate among the models under consideration is obtained as 

% susy (s) + T? usy {s)z + T 2 susy {s)z 2 , 

(3.44) 

where z = cos 9 with 9 being the angle between the momenta of the lepton l + and the 

A b in the center of mass of leptons, v — y 1-is the lepton velocity, A = A(l,r, s) = 

(1 — r — s ) 2 — Ars is the usual triangle function, s = q 2 /m\ b and r = m\/m\ b . The functions 


d 2 T 


SUSY , 


dsdz 


(U«) = 


16384vr 5 


|V«,V;:| 2 tyA(l,r,5) 
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T 0 SUSY (s), T\ SUSY is) and T 2 SUSY (s) are obtained as 

% SUSY (s) = 32m 2 £ m 4 Ab s(l + r - s) ( \V 3 \ 2 + \£ 3 \ 2 ) 

+ 64rrif rn : Ah (1 — r — 3 ) Re V*£ 3 + X> 3 £]* 

+ 64rn Ah \/r'(6rn 2 — m\ b s )Re T>*£\ 

+ 64777,^777,^^12mA 6 sRe V*£ ?J + (1 - r + s)Re X>*X> 3 + £i*£ 3 
+ 32m 2 Ab (2mj + (1 - r + s)m Ab yfr Re + B\B 2 


— tt7a 6 (1 — r — s) Re .4 *£>2 + A\B\ 


— 2 Jr[ Re 


A\B] +m 2 A s Re 


A* 2 B 2 


+ 8m 2 A \4m 2 e (l + r - s) + m 2 Ab (1 - r) 2 - s 2 j ( |A| 2 + \Bi\ 2 ) 

+ 8777a 6 |4777^ A + (1 + r — s)s + m 2 Ab s (1 — rf — s 2 | ( |Al 2 | 2 + \B 2 \ 2 j 
- 8777 ^ 14mf (1 + r-s) -m 2 Ab (1 - rf - s 2 | ( \D 3 \ 2 + |£i| 2 ) 


V* 2 £ 2 


+ 8777| b sn 2 <| — 8m\ b s\/rRe 
- 4(1 - r-s ) Re \v{£ 2 + V\£ x 


+ 4(1 - r + s)v / rRe X>*X> 2 + B\£ 


+ m A , 


(1 -rf-s 2 \V 2 \ z + \£ 2 \ 


- 8m\ b 4t?7£ 


(1 - r) 2 - s(l + r) Re V\K X + £{S X 
+ (4t?7 2 - m\ b s) (1 - r) 2 - s(l + r) ( |d/i| 2 + \Hi | 2 
+ 4m\ b y/rs 2 (4m 2 - m 2 Ab s) Re G 3 H 3 

- 8m Ab s 1 2y/r (4m 2 - m^s) (1 - r + s) Re \G\G 3 + WfU 3 
+ 4m^(l: — r + s)Re T>\K 3 + ^ 3 5 3 + T> 3 JC\ + £ 3 £>i 
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4m*(l - r - s)Re X>JS 3 + ^/C 3 + X> 3 <Si + £*/Ci 
2(1 — r — s) (4m| - m^s) Re C?*% 3 + %*C/ 3 

™a 6 (1 — r) 2 — s(l + r) (|/Ci| 2 + l^l 2 ) j 
32m^ 6 v / ^'§‘|2m^Re T>*S\ + £* X 1C\ + (4 mj — m\ b s) Re Q\T~L\ 

8^a 6 '® 2 ^4\A’ Re K\S\ + 2m\ b y/r(l — r + s)Re /C*/C 3 + <S*<S 3 

2mA 6 (l — r — s)Re /C*<S 3 + <S*/C 3 
(4m^ - mi i s)(l + r - s) ( |£ 3 | 2 + [Rsf ) 


4rrii (1 + r — s)Re 


®S/C 3 + £& 


8mey/rR,e 


+ 8m® s 3 < (1 + r - s) ( |/C 3 | 2 + |S 3 | 2 ) + 4vrRe 


v;s 3 + s:jc 3 


ic;s 3 


(3.45) 


7) 5f/5> (s) = —32m\ b meV\v(l — r)Re [A\Q\ + B\R\ 

-( 


- 16m\ b sv\^l2 R e[AlV 1 )^2 R e(B* 1 £ 1 
+ 2 m Ab R e (B\V 2 - B*V 1 + - A*£ 2 

+ 2m Ab m e R e (a]H 3 + B\Q 3 - A\U X - B* 2 Q X 

+ 32 m\ b s vV^A-f m Ab (l — r)Re(^A* 2 T , 2 ~ B* 2 £- 
+ V^Re (A;V i + A\V 2 - i - 



(3.46) 
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and 


q-SUSY 

• 9 


(s) — —8m\ b v 2 x{^ \Ai\ 2 + \Bi\~ + |"Pi| 2 + |^i| 2 ^ 
+ 8 m ^ b § v 2 X ^ 1 ^ 4 . 21 2 + 1^21 2 + l ^) 2 + |^21 2 ^ • 


(3.47) 


Integrating the Eq.(3.44) over z in the interval [—1,1], we obtain the differential decay 
width only in terms of 3 as 


dTsuSY _ G 2 F a ln m Ab It/ v *\2 /T -r5C/5Y/ S \ , 1 T SUSY (- 
dl ^ ' - 8192;r 5 |V * Vtel ' VA [ T ° (S) + T 2 (S 


(3.48) 


The differential decay rate of RS C is found from dr ^7 sy (g) by replacing C'q,, , Cq 2 and 

Cq 2 with zero. In the case of SM, (5) is found from the supersymmetric differential 
decay rate via setting C ' 7 e , C ' Q e , C[ 0 , Cq 1 , Cq i , Cq 2 and Cq 2 to zero. 


3.2 The differential branching ratio 

In this subsection, we numerically analyze the differential branching ratio that depends on 
q 2 for the A& —> M + £~ decay in SMLCSR, SUSY and RS C scenarios. In order to discuss 
the variation of the differential branching ratio with respect to g 2 , we shall present some 
values of input parameters in table 3 besides the form factors as the main inputs. 



q 2 [GeV 2 /c 4 ] 

Figure 1: The dependence of the differential branching ratio on q 2 for the A & —> 
transition in the SMLCSR and RS C models. The lattice QCD [18] and recent experimental 
data by LHCb [11] Collaboration are also included. 
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Some Input Parameters 

Values 

m tl 

0.10565 GeV 

m T 

1.77682 GeV 

m c 

1.275 ± 0.025 GeV 

m b 

4.18 ± 0.03 GeV 

m t 

173.21 ±0.51 ±0.71 GeV 

mw 

80.385 ± 0.015 GeV 

™a 6 

5.6195 ± 0.0004 GeV 

m A 

1.11568 GeV 

T A b 

(1.451 ±0.013) x 10 - 12 s 

h 

6.582 x 10 ~ 25 GeV s 

g f 

1.166 x 10 " 5 GeV " 2 

^em 

1/137 

IWSI 

0.040 


Table 3: The values of some input parameters used in our calculations, taken generally 
from PDG [42], 



q 1 [GeVV] 


Figure 2: The dependence of the differential branching ratio on q 2 for the A & —y A r + r 
transition in the SMLCSR and RS C models. 

By using all these input parameters and the form factors with their uncertainties, we 
present the dependence of the differential branching ratio of the A& —> Ad + £~ on q 2 in 
SMLCSR, RS C and different SUSY models in figures 1-6. In these figures we also show the 
experimental data provided by LHCb [11] as well as the existing lattice QCD predictions 
[18]. We do not present the results for e in the presentations since the predictions at e 
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q 2 [GeVV] q 2 [GeV 2 /c 4 ] 


Figure 3: The dependence of the differential branching ratio on q 2 for the A& —>■ A/x + /r~ 
transition in SMLCSR and SUSY 1 and If models. The lattice QCD [18] and recent exper¬ 
imental data by LHCb [11] Collaboration are also included. 
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Li 


15 


20 


q 2 [GeV 2 /c 4 ] 


q 2 [GeV 2 /c 4 ] 


Figure 4: The dependence of the differential branching ratio on q 2 for the A& —y A/x + /r _ 
transition in SMLCSR and SUSY 111 and SO(IO) models. The lattice QCD [18] and recent 
experimental data by LHCb [11] Collaboration are also included. 
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q 2 [GeV 2 /c 4 ] q 2 [GeV 2 /c 4 ] 

Figure 5: The dependence of the differential branching ratio on q 2 for the A & —y Ar + r 
transition in SMLCSR and SUSY 1 and 11 models. 
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q 1 [GeV 2 /c 4 ] q 1 [GeV 2 /c 4 ] 


Figure 6: The dependence of the differential branching ratio on q 2 for the A & —y A r + r 
transition in SMLCSR and SUSY Iff and SO(IO) models. 

channel are very close to those of /j. 

From figures 1-6 we see that 

• for all lepton channels, the SMLCSR and RS C models have roughly the same pre¬ 
dictions except for some values of q 2 at which there are small differences between 
predictions of the SMLCSR and RS C models on the differential branching ratio. 

• The areas swept by the SMLCSR are wider compared to those of lattice QCD [18] 
existing in the // channel but they include those predictions. 

• The experimental data in the intervals 4 GeV 2 /c 4 < q 2 < 6 GeV 2 /c 4 and 18 GeV 2 /c 4 
< q 2 < 20 GeV 2 /c 4 cannot be described by the SMLCSR, lattice QCD or RS C mod¬ 
els. In the remaining intervals the SMLCSR, lattice and RS C models reproduce the 
experimental data, except for 6 GeV 2 /c 4 < q 2 < 8 GeV 2 /c 4 , for which the datum 
remains outside of the lattice predictions. 

• In the r channel, the bands of the SMLCSR and RS C scenarios intersect each other, 
except for higher values of q 2 , for which the errors of the form factors do not kill the 
differences between the two model predictions. 

• At all lepton channels, the SUSY models show overall considerable deviations from the 
SMLCSR, lattice QCD and experimental data although they include the predictions of 
these models for some values of q 2 . The maximum deviations of the SUSY predictions 
from the results of SMLCSR, lattice QCD and experiment belongs to the SUSY II 
such that the SMLCSR, lattice QCD and experimental results remain out of the 
regions swept by the SUSY II model at higher values of q 2 . 
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• In the n channel, the experimental data in the interval 18 GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 
are reproduced by SUSY I, III and SO(IO) but not by SUSY II. Note that in this 
interval other models (SMLCSR, lattice QCD and RS C ) also can not describe the 
experimental data. 

• Again in the /i channel, the experimental data in the interval 4 GeV 2 /c 4 < q 2 < 6 
GeV 2 /c 4 cannot be reproduced by any SUSY models like the SMLCSR, lattice QCD 
and RS C scenarios. 

• In the case of r as the final lepton, there are considerable differences between different 
SUSY models’ predictions and that of the SMLCSR and these cannot be completely 
killed by the errors of form factors. The maximum deviations of the SLISY results 
from the SMLCSR predictions belong to the SUSY II at higher q 2 values. 


3.3 The lepton forward-backward asymmetry 


In this subsection, we would like to present the results of the lepton forward-backward 
asymmetry obtained in different scenarios. The lepton Afb is defined as 


Afr(s) — 


I' <* 2 r , J f° cPr , J 

/ -TT7:(2,f>)<I2 - / -FTJlC.SJ'fe 
I o dsdz i cLsdz 

f 1 d 2 T . ^ , f° d 2 T . ^ , 

/ yrr Z { z , s )dz+ / —— (z,s)dz 

/q dsdz J_ i dsdz 


(3.49) 



q- [GeV 2 /c 4 ] 


Figure 7: The dependence of the Afb on q 2 for A& —> A/i + yU transition in SMLCSR, lattice 
QCD [18] and RS C models together with recent experimental data by LHCb [11]. 
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q 1 [GeVV] 


Figure 8: The dependence of the Afb on q 2 for A & —* A r + r transition in SMLCSR and 
RS C models. 
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Figure 9: The dependence of the Afb on Q 2 for A-b A/i + /x transition in SMLCSR, lattice 
QCD [18] and SUSY I and II together with recent experimental data by LHCb [11]. 
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Figure 10: The dependence of the Afb on q 2 for A & —> A/d + n transition in SMLCSR, 
lattice QCD [18] and SUSY 111 and SO(10) together with recent experimental data by 
LHCb [11], 
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Figure 11: The dependence of the Afb on q 2 for A& —>■ Ar + r transition in SMLCSR and 
SUSY I and If scenarios. 
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Figure 12: The dependence of the Afb on q 2 for A b —* Ar + r transition in SMLCSR and 
SUSY Ill and SO(10) scenarios. 


Considering the form factors with their uncertainties from [19], we plot the dependence 
of the lepton forward-backward asymmetry on q 2 for the decay under consideration in both 
lepton channels in the SMLCSR, RS C and different SUSY models in figures 7-12. From 
these figures, we obtain that 

• in the fi channel, the SMLCSR, lattice QCD and RS C models predictions on Afb 
coincide with each other. Except for the lattice QCD, they can only describe the 
experimental data existing in the 0 GeV 2 /c 4 < q 2 < 2 GeV 2 /c 4 and 18 GeV 2 /c 4 
< q 2 < 20 GeV 2 /c 4 regions. The remaining data lie out of the swept regions by all 
these models. 

• As far as the SLISY models are considered, in the // channel, the SUSY models have 
predictions that deviate from the SMLCSR and lattice QCD predictions, considerably. 
All SUSY models reproduce the experimental data in the regions 0 GeV 2 /c 4 < q 2 < 2 
GeV 2 /c 4 and 18 GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 . The other data also remain out of the 
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regions swept by different SUSY models except for the SUSY II, which reproduces 
the experimental data also in the interval 15 GeV 2 /c 4 < q 2 < 18 GeV 2 /c 4 . 

• In the case of the r lepton, the SMLCSR and RS C have roughly the same predictions 
on Afb- 

• In the r lepton channel, the SMLCSR and SUSY I have roughly the same predictions 
for Afb ; however, the remaining SUSY models’ predictions deviate from the SMLCSR 
predictions considerably, although they intersect each other at some points. 


4 Conclusion 

In the present work, we have analyzed the scmileptonic A h —> Ai. + t~ decay mode in 
SMLCSR, different SUSY models and the RS C scenario. Using the form factors calculated 
in light cone QCD sum rules in the full theory [19], we evaluated the differential branching 
ratio and lepton forward-backward asymmetry for different leptons in those scenarios. We 
also compared the results obtained via SMLCSR, RS C and different SUSY scenarios with 
the recent experimental data provided by LHCb [11] as well as the existing lattice QCD 
predictions [18] on the considered quantities. We observed that the regions swept by the 
SMLCSR model include the RS C predictions although they are somewhat wider compared 
to those of RS C models for the considered physical quantities. The SMLCSR predictions on 
the considered quantities in the present work are overall consistent with the lattice QCD 
predictions provided by Ref. [18]. 

The predictions of different SUSY models on the differential branching ratio deviate 
considerably from the SMLCSR and lattice predictions. The maximum deviations belong 
to the SUSY II model. In the case of Afb and the /j channel, the predictions of different 
SLISY models have considerable deviations from the SMLCSR and lattice QCD predictions. 
For Afb and the r channel, the SLISY I and SMLCSR have roughly the same predictions 
but the other SUSY models have predictions different from that of the SMLCSR. 

The experimental data on the differential branching ratio in the /i channel can be re¬ 
produced by SMLCSR, lattice QCD and RS C models except for the intervals 4 GeV 2 /c 4 

< q 2 < 6 GeV 2 /c 4 and 18 GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 , which cannot be described by 
SMLCSR, lattice QCD or RS C models. As far as the SUSY models are considered, differ¬ 
ent SUSY models also cannot reproduce the experimental data in the interval 4 GeV 2 /c 4 

< q 2 < 6 GeV 2 /c 4 . However, except for SUSY II, the remaining SUSY scenarios can explain 
the experimental data in the region 18 GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 . 
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In the case of Afb and the /i channel, the SMLCSR, RS C and different SUSY models 
can only describe the experimental data existing in the 0 GeV 2 /c 4 < q 2 < 2 GeV 2 /c 4 and 
18 GeV 2 /c 4 < q 2 < 20 GeV 2 /c 4 regions. The other existing data remain out of the swept 
areas by these models, except for SUSY II, which can also reproduce the experimental data 
in 15 GeV 2 /c 4 < q 2 < 18 GeV 2 /c 4 . 

More experimental data in the /i channel related to different physical quantities asso¬ 
ciated with the At —» A p + mode, the future experimental data in the r channel and 
comparison of the results with our predictions on the quantities considered in the present 
work may help us in the course of searching for new physics effects. 
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